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We theoretically study how to control transport, bound states, and resonant states of a single photon in a 
one-dimensional coupled-cavity array. We find that the transport of a single photon in the cavity array can be 
controlled by tuning the frequency of either one or two cavities. If one of the cavities in the array has a tunable 
frequency, and its frequency is tuned to be larger (or smaller) than those of other cavities, then there is a photon 
bound state above (or below) the energy band of the coupled-cavity array. However, if two cavities in the array 
have tunable frequencies, then there exist both bound states and resonant states. If the frequencies of the two 
cavities are chosen to be much larger than those of other cavities, and the hopping couplings between any two 
nearest-neighbor cavities are weak, then a single photon with a resonant wave vector can be trapped in the region 
between the two frequency-tunable cavities. In this case, a quantum supercavity can be formed by these two 
frequency-tunable cavities. We also study how to apply this photon transport control to an array of coupled 
superconducting transmission line resonators. 

PACS numbers: 03.67.Hk, 03.65.-w, 05.60.Gg 



I. INTRODUCTION 

In a quantum network based on photons JTj], the nodes can 
be regarded as information processing stations (e.g., quantum 
computers), while the links between any two nodes are pro- 
vided by the information carriers (e.g., photons). Due to the 
high-speed transmission and low dissipation in optical fibers, 
photons are considered to be excellent information carriers 
(both for classical and quantum transmissions). 

In recent years, with the development of nano-optics, 
numerous photonic information processing proposals have 
been implemented by using on-chip solid state devices, such 
as semiconducting microcavities (e.g., Ref. [2]) and super- 
conducting transmission line resonators (e.g., Refs. Dj-Jal). 
Therefore, how to realize on-chip single-photon devices (e.g., 
Refs. 1I9 UI1I0 becomes now an increasingly important re- 
search area. For example, single-photon switches (e.g., 
Refs. (l2|4ldj), which control single-photon transport at will 
(e.g., Refs. Ill7l [l8ll ). play an important role in this area. 

Coupled-cavity arrays (CCAs) II 19143011 are one type of pho- 
tonic system, which has been proposed to process photonic 
quantum information. Compared with the usual waveguides, 
which have only a linear dispersion relation, the nonlinear dis- 
persion of CCAs can result in the emergence of bound states 
of single photons. Many proposals have been put forward to 
realize quantum switches in CCAs, which could be used to 
control single-photon transport. For example, the reflection 
and transmission of single photons in a one-dimensional cou- 
pled resonator waveguide can be controlled by a tunable two- 
level system inside one of the cavities [14]. Moreover, con- 
trollable single-photon transport in a one-dimensional CCA 
with a tunable hopping coupling has recently been stud- 
ied £01. 



In this paper, we study another approach to realize con- 
trollable single-photon transport in a one-dimensional CCA, 
considering either one frequency-tunable cavity (FTC) or two 
FTCs. This work is motivated by recent mostly experimen- 
tal results on frequency-tunable transmission line resonators 
(e.g., Refs. Il3ll437ll ). where the frequencies of the resonators 
can be changed by varying either the boundary condition of 
the electromagnetic wave or the magnetic flux through the 
SQUIDs used to construct the transmission line resonators. 
In contrast to Refs. lfl4H27ll . here the photon transport is con- 
trolled by tuning the frequency of the cavity, and there is no 
additional two-level system, placed inside one of the cavities, 
to control photon transport. Therefore, this proposal seems to 
be simpler and easier to implement experimentally than those 
in Refs. iTfl l27ll . By changing the frequency of either one 
FTC or two FTCs, the reflection and transmission of a single 
photon in the coupled-cavity array can be controlled. We also 
study the photon bound states and photon resonant states Ik38ll 
in this coupled-cavity array. 

For the coupled-cavity array with one frequency-tunable 
cavity, if the frequency of the frequency-tunable cavity is 
larger than that of other cavities, there is a bound state above 
the energy band of the corresponding bosonic tight-binding 
model [39]; while the bound state is below the energy band 
when the frequency of the frequency-tunable cavity is smaller 
than that of other cavities. 

For the CCA with two FTCs, we find that there exist bound 
states around the FTCs. Moreover, when the frequencies of 
the two FTCs are much larger than those of other cavities and 
the hopping couplings between any two nearest-neighbor cav- 
ities are weak, for resonant wave vectors, a single photon can 
be in resonance with the CCA and then remain trapped in the 
cavities between the two FTCs. A single photon in resonance 



with the CCA behaves as a photon inside a supercavity |27 

This paper is organized as follows: In Sec.[TTJ we study con- 
trollable single-photon transport and bound states in a CCA 
with one FTC. In Sec. [TTTJ we study controllable single-photon 
transport, bound states, and resonant states in the CCA with 
two FTCs. In Sec. [IV] we present a possible experimental 
implementation of our proposal using superconducting trans- 
mission line resonators. A summary is given in Sec.lVl 



H. COUPLED CAVITY ARRAY WITH ONE 
FREQUENCY-TUNABLE CAVITY 
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As schematically shown in Fig. [T] we consider a one- 
dimensional coupled-cavity array, which consists of a chain 
of N cavities. We assume that N is a large enough number, so 
periodic boundary conditions become reasonable. For speci- 
ficity, and without loss of generality, we assume that N is an 
odd number. The distance between any two nearest-neighbor 
cavities is do- The cavities, except the central one (i.e., the 0th 
cavity), have the same resonant frequency oj c . The central 0th 
cavity has the resonant frequency (1 + A)o> c , where A is used 
to characterize the detuning between the 0th cavity and other 
identical cavities, and assume to be varied for controlling the 
photon transport properties in this system. The frequency of 
the 0th cavity can be larger (A > 0) or smaller {A < 0) than 
those of other cavities. Any two nearest-neighbor cavities are 
coupled via a homogeneous hopping interaction of strength /. 
The Hamiltonian (with ft — 1) of the CCA reads 

A-l AM 

Hqq A - Acj c a^ao + ^ Lo c a.aj- ^ J{a\cij+\ + a. +l cij), 

:_ JV-1 N-\ 
J~ 2 ' 1 

(!) 

where aj (at) is the annihilation (creation) operator of the jth 

cavity. The superscript "1" in 6qqa denotes that the CCA 
contains one frequency-tunable cavity. The first two terms 
in Eq. (HJ are the "free Hamiltonian" of the CCA, while the 
last term in Eq. ([D represents the hopping interaction, with 
strength J, between any two nearest-neighbor cavities. For in- 
stance, the term aAj+x means that a photon is annihilated in 
the ( j + 1 )th cavity and another photon is created in the jth 
cavity. Hereafter, we only use instead of the sum shown 
in Eq. (fTJ. Since the frequency, (1 + A)cl> c , of the 0th cavity 
should be nonnegative, then A > — 1. For A — 0, the above 
Hamiltonian ([TJ reduces to the usual bosonic tight-binding 
(btb) Hamiltonian 



Hbtb = oj c V fljfl j - J ^Ja\a 



aj), 



which can be diagonalized 

#btb = Qk<i\ak 
k 

by using the discrete Fourier transform 

a k = — p= y. exp(ikjdo) a. 
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FIG. 1: (Color online) Schematic diagram of a one-dimensional 
coupled-cavity array (CCA): (a) a lattice model for the CCA. Each 
disk represents a cavity (labeled by the integer below each dot) with 
frequency shown right above each cavity. The larger disk (located at 
the site j = 0) represents the frequency-tunable cavity; (b) the energy 
levels of the CCA; (c) coupled superconducting transmission line res- 
onator array. The central cavity located at the site j = is shown in 
orange. This central cavity has a tunable frequency (1 + A)ll> c , where 
A is the detuning parameter used to control the transport of single 
photons through this coupled-cavity array. If A > (A < 0), then 
there is a photon bound state above (below) the energy band of the 
coupled-cavity array. Thus, the central cavity acts as a tunable "im- 
purity". The arrows on the far left (right) schematically indicate the 
incoming (outgoing) photons. The blue dots represent the remaining 
cavities which are not shown here. The integers below each cavity 
label each one of them. 



and the periodic boundary conditions, where 
Elk — oj c - 27 cos(faafo) 



(5) 



is a nonlinear dispersion relation, which is an energy band. 
Hereafter, the distance do between two neighboring cavities is 
scaled as unity, and k - Innu/N (with —N/2 < n^ < N/2) are 
the photon wave vectors. 



A. Controllable single-photon transport for one 
frequency-tunable cavity 

Since the total excitation number operator = Yjj a )&j °f 

the CCA is a conserved observable, i.e., [N, H^.q A ] = 0, it 
is reasonable to restrict our discussions to the single-particle 
excitation subspace for studying single-photon transport. A 
general state in the single-excitation subspace can be written 

as 



M = ^>|1 ; .>, 



(6) 
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where the state \ = |0) ® ■ ■ ■ ® |1); ® ■ ■ ■ ® |0) rep 
the case when the jth cavity has one photon, while oth 
ities have no photons. Also, Cj is the probability ampli 
the state Using the discrete scattering method stu 
Ref. 111411 and the eigenequation dilAoj) = lo\(l>}, we c 

- J{cj+\ - Cj-\) = [a>- (1 + A5jo)u> c ]cj, 

where Sjo is the Kronecker delta function. 

Without loss of generality, we assume that a single 
with frequency a> — £2^ is injected from the left side 
CCA, and then the photon probability amplitudes c> 
sumed to have the following solutions, 



j<0, 
j>0, 



where r and s are the photon reflection and transmissi 
plitudes, respectively. Here "i" denotes the imaginary u 
cept when specified otherwise. It is easy to check that 
is the solution of Eq. (O when j + 0. Connecting Eq. i|7]i ai 
j — with the continuity condition, 1 + r = s, for the wave 
function, we obtain the photon reflection amplitude 



Aa> c 



2iJ sin k - Acl> c ' 
which leads to the photon reflection coefficient 

(Aco c ) 2 



R{k,A) = \r\- 



4J 2 sin 2 k + {Acj c ) 2 ' 



(9) 



(10) 



The reflection coefficient has three symmetric relations: 
R(k,A) = R(-k,A), R(tt/2 - k,A) = R(n/2 + k,A), and 
R(k, A) = R(k, -A). In Fig.|2] the reflection coefficient R(k, A), 
as a function of the detuning parameter A, is plotted for 
k = 0.01, 7r/8,7r/4, and n/2. It can be seen from Fig. |2]that the 
photon reflection coefficient R(k, A) can be tuned from zero to 
one by changing the detuning parameter A. 



B. Bound states for one frequency-tunable cavity 

A bound state may be formed when a particle is scattered 
by a localized potential. In the present model, the on-site extra 
energy Aa> c acts as a potential. Therefore, a priori, this system 
may have bound states. First, we give a qualitative analysis 
of the bound states in the CCA. We now apply the discrete 
Fourier transform, defined in Eq. ©, to express the Hamilto- 
nian (Q]l in wave vector space as follows, 
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In terms of the Hamiltonian in Eq. ( fTTT i. we obtain the Heisen- 
berg's equation of motion for the operator a^ 



ia^ — ^l^a^ 



Aa> c 
~~N~ 



Z 



(12) 
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FIG. 2: (Color online) The photon reflection coefficient R{k, A) ver- 
sus the detuning parameter A of the Oth cavity is plotted for k = 
0.01,^/8,^/4 and n/2, where the parameters are taken in units of 
co,-, and J/co c = 0.01. Recall that the on-site detuning of the central 
(Oth) frequency-tunable cavity is Aa> c . A relatively small amount of 
detuning can make the reflection coefficient near one. 



For small A, we can assume that ia^ = oj^ak- By introducing 
the operator b = 2a Sjt, we obtain 



A(jJ r 



1 



a k = 



N a>k-Qk 



Making the summation 



L-i m Z_ i 



N ti>k-Qk 



then the frequencies Uk are determined by the equation 



k 



l 



cj k - D.k 



= 1. 



(13) 



(14) 



(15) 



Equation Sl5[ can be solved numerically. In Figs. [3] the func- 
tions f\(to) = Acl) c 2jt[l/(w-Qjt)]//V and/2(w) = 1 are plotted 
for A > and A < 0. The values of a> corresponding to the 
crossing points of both curves f\(cS) and fiito) are the solu- 
tions cok that satisfy Eq. d!51 l. Obviously, when A > there is 
a bound state above the energy band, while for A < there is 
a bound state below the energy band. These bound states are 
shown as black circles in Fig. [3] 

Below, we analytically study the bound states. Since we 
choose the frequency-tunable cavity as the coordinate origin, 
the Hamiltonian (Q~|) of the system is symmetric around the 
Oth cavity. The eigenstates of the Hamiltonian (HJ have ei- 
ther symmetric or asymmetric parities. For the asymmetric 
case, we have the relation c ; = -C-j, which implies cq = 0. 
Therefore, the frequency change of the frequency-tunable cav- 
ity will not affect the asymmetric eigenstates. For the symme- 
try case, we assume the following solution 



Cj = A\n\K 



(16) 
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FIG. 3: (Color online) The functions fi(a>) = Auj c 2*[l/(<^ - Q t )]/N 
(blue) and fi(oS) = 1 (red) versus the scaled frequency (j/co c . The 
values of a) corresponding to the crossing points between the curves 
for the functions fi{oS) and foicD) are the frequencies 0)k- In other 
words, the oi^'s satisfy f\(b)k) = fzU^k)- Here, the parameters are 
chosen as N = 21, 0) c = 1, Jju) c = 0.01. Note that (a) and (b) use 
A = 0.2 and A = -0.2, respectively, and this is their only difference. 
Obviously, for A = 0.2 (A = -0.2), there is a bound state (shown as a 
black circle) above (below) the energy band (0.98 < co/a) c < 1.02). 



where n is a parameter introduced to describe the bound state 
of the Hamiltonian (|T). 

Substituting the solution ([Tol l into Eq. ©, we obtain 

- Jfi 2 + Acd c {i + J = 0. (17) 

Equation (jTTJ has two solutions 

_ -Aco c + V47 2 + (Acj c ) 2 
ix ± - — . (18) 

When A > 0, we choose the solution while for the case 
A < we choose the solution yu_. The corresponding eigen- 
frequencies are 

oj ± = cj c ± y/4J 2 + (Ato c ) 2 . (19) 

The relations <d+ > co c + 2J and <y_ < u> c - 2J mean that 
the bound state is above and below the energy band, respec- 
tively, as shown in Figs. [3|a) and Ob). These analytical re- 
sults are consistent with our previous analysis in Eqs. dTTT - FT3T > 
and Fig. [3] Note that /z+(|/l|) = - fJ~(— \A\), so |c,-| 2 is the same 
for the two bound states. In Fig.|4j we plot |c,-| 2 as a function 
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FIG. 4: (Color online) Photon probability |c/| 2 versus the location j 
of the cavity. The c/s are introduced in Eq. ©. Other parameters 
are taken as a> c = 1, \A\ = 0.2, and J/w c = 0.01. In this case, since 
\cj\ 2 is peaked at the cavity located at j = 0, then a single photon is 
localized around the central (0th) frequency-tunable cavity. 

of the lattice parameter j. Figure H] shows that a single photon 
is mainly localized around the central (0th) frequency-tunable 
cavity. 

C. Links to localized excitations in solids 

Periodic solid state systems exhibit bands. Adding local- 
ized defects to a translationally invariant structure induces lo- 
calized states around those defects. In general, adding a defect 
(i.e., anything that breaks translational symmetry, like an im- 
purity, or an interface) is enough to create gap states outside 
the continuous bands. Thus, gap states outside the bands are 
linked to localized states. 

Figure [3] in this work is related to Fig. 9.9 in page 395 of 
Madelung's classic textbook J40ll on solid state physics. In- 
deed, equation (9.37) in ll40Tl is related to Eq. ( fT3] l in our pa- 
per. The links between them is that a defect added to a peri- 
odic structure tends to localize excitations around the defect, 
and this localized state corresponds to gap states. 



III. COUPLED CAVITY ARRAY WITH TWO 
FREQUENCY-TUNABLE CAVITIES 

As schematically shown in Fig. [5] we now consider the case 
when there are two FTCs in the CCA. These two FTCs are 
located in the -fifth and the dth cavities, respectively. The 
Hamiltonian can now be written as 

^cca = ^iCL> c a[ d a-d + A2io c a ] d ad + co c ^ a]aj 

j 

-J^(a)a j+l +a\ +1 aj), (20) 
j 

where Ai and A2 are, respectively, used to describe the fre- 
quency of the -dth and the dth cavities, and J is the hopping 
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FIG. 5: (Color online) Schematic configuration of a one-dimensional 
coupled-cavity array with two frequency- tunable cavities: (a) a lat- 
tice model for the CCA. Each disk represents a cavity (labeled by 
the integer below each cavity) and with frequency shown right above 
it. The two larger disks (located at the sites j = ±d) represent the 
two frequency-tunable cavities; (b) the energy levels of the CCA; (c) 
schematic diagram of a coupled superconducting transmission line 
resonator array. The two detuned cavities have frequencies (1 + A{)u) c 
and (1 + A2)(x) c . The two detuning parameters, Ai and Ai, control the 
transport of single photons in the array of coupled cavities. The in- 
coming photon can be localized at the cavities located at j = ±d, 
forming localized bound states (in Fig. [8). Alternatively, when the 
photon wave vector is in resonance, the photon can be confined in 
the region between these two cavities, as shown in Fig. [9] 
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FIG. 6: (Color online) Photon reflection coefficient R plotted versus 
the detuning parameter A for k = O.ln, 0.2n, and 0.37T, where co c = 
1, d = 5, and J/uj c = 0.01. Here we consider two frequency-tunable 
cavities located at j = —d and j - d. A relatively small amount of 
detuning can make the photon reflection coefficient near one. 



Eq. (l2TT i is assumed as 



[ e ikj + re-' kj , 
cj = I Ae lk > + Be'^j, 
\ se ik K 



j < ~d, 
-d < j < d, 
j > d. 



(22) 



Substituting the solution ( l22b into Eq. ( |2TI ) and using the con- 
tinuity condition at j = -d and j — d 



e -ikd + re ikd = Ae -ikd + Be ikd 



Ae m + Be~ iU 



Jkd 



(23a) 
(23b) 



coupling strength between any two nearest-neighbor cavities. 

~ (21 

The superscript "2" in H C {, A means that the CCA contains two 
frequency-tunable cavities. 



A. Controllable single-photon transport for two 
frequency-tunable cavities 



For the CCA with two frequency-tunable cavities, the total 
excitation number operator = £j a\<Xj is also a conserved 
observable. Similar to Eq. ©, the eigenstates of the system 
can now be written as \a>) = Y*j c Mj)> where the coefficients 
Cj are determined by the equation 

- J(Cj-\ + cj+i) = \cj - co c (l + MS-dj + MSd.j)] cj. (21) 



we obtain the photon reflection amplitude 

r = A Q u c e-^ k + e 4ikd ) (e 2ik - l) 
+ Ao<o c (e 4ikd -l)e ik ] 
x {j (e 2ik (e 2ik - l) - 2A aJ c e ik ] 



„2ik\ 



(24) 



H 2 is 



In Fig. |6l the photon reflection coefficient R(k, Aq) 
plotted as a function of the parameter Aq for different wave 
vectors k = 0.17T, Q.2n, Q.3tt. Figure|6]shows that the reflection 
coefficient R(k, Aq) can be tuned from zero to one by changing 
the detuning parameter Aq. 



B. Bound states and resonant states for two frequency-tunable 
cavities 



For simplicity, and without loss of generality, below we as- 
sume that the parameters for the two frequency-tunable cavi- 
ties are identical, i.e., A[ = A2 = Aq. 

For a photon with frequency co = Q.^, the functions e\p(ik j) 
and exp(-ikj) are the solutions of the Eq. ( l2"TT i when j + -d 
and j + d. Therefore, the general form of the solution for 



To study the bound states of the CCA with two frequency- 
tunable cavities, we assume 



= I Be ik i + Ce~ ik j, 
De ik i, 



j < ~d, 
-d < j < d, 
j>d. 



(25) 
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From Eq. (|2TT > at j = -d, and the condition for the continuity 
of the wave function at j = -d, (namely: Ae 



,kd 



Be 



-ikd 



Ce ), we obtain 



i2J sin k 



B Aqlo l 

C (i2J sin k - Aqu> c ) 
~B ~~ 



AqCO, 



exp(-2ikd), 

exp(-2ikd). 



(26a) 
(26b) 



From Eq. (l2TT i at j = d, and the continuity condition at j - d, 
(that is Be m + Ce~ w = De ikd ), we obtain 



AqO) c 



C _ 

B i2J sin k - AqoJc 

D -i2J sin k 

B i2J sin k - Aqu> c 



exp(2ikd), 



(27a) 
(27b) 



Obviously, the right hand sides of Eqs. d26bb and d27ab should 
be equal; then we obtain 

(Aqco c ) 2 exp(4ikd) = (Aqu> c - i2J sin A:) 2 , (28) 



which implies 



exp(2ikd) 



ill 

Aoo> c 



■ sink 



(29) 



This Eq. d29l is important here because it determines the wave 
vector of the photon states (either bound states or resonant 
states). Below, we will discuss the existence of bound states 
in this system for the two cases shown in Eq. ( 1291 ). These will 
be denoted as positive root and negative root, respectively, of 
Eq. d28| i. Also note that Aq takes two values (±|/lo|) because 
the frequency detuning of the cavities at j = +d can be either 
positive or negative. 

1. Positive root of Eg. \29\ 

We now assume that the wave vector k — x + iy, with y > 0. 
It should be pointed out that the parameter y can take the value 
y — to include the possibility for the existence of some in- 
teresting resonant states. Since the energy of the photon must 
be a real parameter, (i.e., cos k = cos jt cosh y - / sin x sinh y 
is real), this condition can be satisfied in the following three 
cases: 

(1) x = 2nn, with n e Z (hereafter, Z denotes the set of 
integers). In this case, y is determined by the equation 



exp(-2c/y) = 1 + 



2J 

Aqco c 



sinhy. 



(30) 



The solutions of the transcendent equation (l30l > (and other 
transcendent equations below) are determined through the nu- 
merical method briefly sketched in Fig. [7] 

The coefficient relations in Eqs. d26l i and ( l27l i are 



j = ~jj = exp(2dy)- 



1, - 



-1. 



(31) 



When Ao > 0, from Fig. |2ta), we can see that Eq. ( f30b has 
only a "zero solution", y — 0, then A = D — and C - -B, so 
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FIG. 7: (Color online) The functions exp(±2dy) and 1 ± 
2J sinh y/(\Ao\aj c ) are plotted as a function of y. The wave vector 
k = x + iy, with y > here. The intersection points shown inside the 
black circles give the nonzero solutions of the transcendental equa- 
tion \29\ . The parameters used here are a> c = 1, J/a> c = 0.01, d = 5, 
and \Aq\ = 0.2. Equation \29\ determines the wave vector of the pho- 
ton states (either bound states or resonant states). 



the corresponding wave function c ; = 0. When Aq < 0, from 
Fig.|7tb), we know that Eq. (1301 has two solutions: one is zero, 
and the other one is a positive number denoted by yo > 0. For 
the zero solution, the wave function is cj = 0. For the solution 
>'o > 0, we have A = -D = (e 2dy ° - V)B and C = -B. The 
corresponding wave function ( |25| ) becomes 



Ci = 



B{e ld y» - V)e yaj ., 
-2B sinh(y j"), 
-B(e 2d >° - l)e 



yoj 



j < -d, 
-d < j < d, 
j > d, 



(32) 



where B is determined by the normalization condition. Obvi- 
ously, the wave function ( l32l is asymmetric, i.e., Cj - —C-i, 
so it is an odd-parity state. Using the parameters oj c = 1, 
J/a) c = 0.01, d = 5, and \A \ = 0.2, we obtain y = 2.998. The 
photon probability corresponding to the wave function (T32b 
is plotted in Fig. [8] which shows that the probability to find a 
single photon d32l is mostly around the two frequency-tunable 
cavities. 

(2) x = (2n + l)7r, with n e Z. In this case, y is determined 
by the equation 



exp(-2afy) = 1 



2J 

AfjUc 



sinhy, 



(33) 



and the coefficient relations in Eqs. d26l i and ( 1271 ) are the same 
as those in Eq. ( |3TT >. When Aq < 0, from Fig.|7Ja), we can see 
that Eq. (1331 has only a zero solution, y = 0, then A — D — 
and C = -B. The corresponding wave function is cj = 0. 
When Ao > 0, from Fig.|7tb), it can be seen that Eq. (|33T > has 
two solutions, one is zero, and the other is a positive number 
denoted by yo > 0. For the zero solution, the wave function is 
cj = 0. For the solution y > 0, then A = -D = (e 2dyo - l)B 
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FIG. 8: (Color online) Photon probability \cj\ 2 of each coef 
of the normalized wave function given by Eqs. ((6) and ( 132b 
parameters chosen here are (j c = 1, J/to c = 0.01, d = 5, y = 
and with a detuning \Aq\ = 0.2 for the cavity at j = ±d. The f 
is now localized and bound at the two detuned cavities Iocs 
j = ±d = ±5. 



and C - —B. The corresponding wave function become; 



B(e ldya - Y)e i{ln+X)n, e y<s i , 

-2Be i(2n+1) *i swh(yof), 
-B(e 2dy ° - l)e ! '( 2n+1 > r ^e _ >o7 ) 



—d < j < d, 
j > d. 



(34) 



The wave function in Eq. fl34"l i is asymmetric for /', so it is 
an odd-parity state. The square of the module of the wave 
function in Eq. d34"t is the same as that of the wave function in 
Eq. (H2>. 

(3) y = 0. The motivation for studying the case y = 
is to investigate whether there exist some interesting resonant 
states. In this case, x is determined by the equation 



cos(2dx) = 1, 
■ sin x — 0. 



AqCO c 



(35a) 
(35b) 



The solution of Eq. J35at is x — mn/d, with m e Z. For 
Eq. ( I35bl ). the solutions are x = In, with / e Z, for a general 
2J/(Aqoj c ), or any jc when Aqcj c » 27, i.e., 2J/(Aou> c ) ~ 0. 
Connecting the solutions for the two Eqs. (I35ab and ( I35bb . we 
obtain two solutions for the case of y — 0: (a) x = for, then the 
coefficient relations in Eqs. (|26| | and d27b become 



5 



~B 



0, 



c 

5 



-1. 



(36) 



Therefore, the wave function is c, = 0; (b) when Aqu » 27, in 
this case, we choose x = mn/d, then the coefficient relations 
in Eqs. (|26| | and (l27l i are the same as those in Eq. (T36b . The 
wave function becomes 



0, 

2«Z?sin(=), 
0, 



j < ~d, 
-d < j < d, 
j > d, 



(37) 



Obviously, the resonant state (|37] i is an odd-parity state with 
parameter j. The square of the module of the wave func- 
tion in Eq. ( f37T > is plotted in Fig. |9ja). We find that, when 
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FIG. 9: (Color online) Photon probability |c,| 2 , given in Eqs. Q7t 
and J45b . versus the location 7 of each cavity are plotted in figures 
(a) and (b), respectively. Here the parameters are set as co c = 1, 
m = 2, d = 5, and 2J/(A u> c ) ~ 0. Now the photon is not bound at 
j = +d, as in Fig. [8] but it is confined in between j = ±d = ±5. This 
is because a single photon with wave vector (a) 2n/5 and (b) 7r/2 is 
in resonance with the coupled-cavity array between the two detuned 
cavities located at j = ±d = ±5. These are resonant states. 



Aqoj c » 27, (i.e., the frequencies of the -dth and dth cavities 
are very largely detuned from those of other cavities and the 
hopping coupling 7 between two nearest-neighbor cavities is 
weak), the photon with wave vector mn/d (m e Z) can pro- 
duce a resonance in the region between the two FTCs, once it 
is injected there. Therefore, this resonant photon state is con- 
fined between the two cavities located at j = +d, as shown in 
Fig. Sa). 



2. Negative root of Eq. \29\ 

We now write the wave vector as k — x + iy, with y > 0. 
Since the energy of the photon must be real, (i.e., cosfc = 
cos x cosh y - i sin x sinh y should be real), this condition can 
be satisfied in the following three cases: 

(1) x = 2nn. In this case, y is determined by the equation 



- exp(-2dy) = 1 + 



27 

Aq(a) c 



sinhy, 



(38) 



and the coefficient relations in Eqs. d26i i and ( f2Tb become 



A _ D _ 
B ~ ~B ~ e 



2dx 



+ i, £ = 1. 

B 



(39) 



In this case, when Aq > 0, we know from Fig. [T^c) that 
Eq. ( f3~8b has no positive solution. When Aq < 0, from 



Fig. |3d), we know that Eq. (1381 has a positive solution y\ 



Then A - D - (e zdV| + 1)5 and C = B. Therefore, the w 
function becomes 



B(e 2dy > + l)e»j, 
2B cosh(viy'), 

B(e 2d » + l)e~ yij 



j < ~d, 
—d< j < d, 
j > d. 



The wave function in Eq. d40l > is an even-parity state. U 
the parameters oj c = 1, J/oj c = 0.01, d = 5, and \Aq\ = 0.2. 
obtain yi = 2.998. We found that the difference betweei 
and y\ is very small, on the order of 10~ 12 . Thus the ph( 
probability |c,-| 2 here looks like the one in Fig. [8] 



(2) x = (2n 
equation 



+ 1)tt. In this case, y is determined by 



- exp(-2dy) = 1 - 



27 

/low f 



sinhy, 



and the coefficient relations in Eqs. d26l i and ( T2"7b are the s; 
as those in Eq. ( |39l . When Ao < 0, from Fig.|7Jc), we can 
that the Eq. (PiTT i has no positive solution. When An > 0, f 
Fig- Eld), we know that Eq. ( |4T1 i has one positive solutior 
Then A - D - (e 2dyi + l)B and C = B. So the corresponc 
wave function now becomes 



2Bep"+V*i coshCyij), 
^(e 2 ^ + l)e ,(2 " +1)7r >e- 



7 < -</, 
-rf < j < d, 
j>d, 



This wave function (l42l i is an even-parity state. The squar 
the module of the wave function in Eq. (l42l is the same as 
of the wave function in Eq. ( f40b . 

(3) y = 0. In this case, x is determined by the equation 



cos(2c/x) = — 1, 
■ sinx = 0. 



Aq(lI c 



(43a) 
(43b) 



The solutions of Eq. ( 143 at are x = (2m + l)n/(2d), and the 
solutions of Eq. ( I43bb are x = or any x when /Iqw c » 
27. Therefore, the solutions meeting the two equations ( I43ab 
and (I43bb at the same time are x = (2m + l)n/(2d), when 
Ai)Ci) c » 27. Then the coefficient relations in Eqs. d26l ) and 
(ITTT i become 



AD C 

- = - = 0, - = 1. 

B B B 

Therefore, the wave function now becomes 



0, 



2Bcos(i 2 ^M), 
0, 



j < ~d, 
—d < j < d, 
j > d. 



(44) 



(45) 



This resonant state d45b is an even-parity state. The square of 
the module of the wave function d45b is plotted in Fig. [9fb). 
This figure shows that a single photon with wave vectors 
(2m + l)n/(2d) can be in resonance in the region between the 
two frequency-tunable cavities. For the two cases of resonant 
states given in Eqs. (|37| | and (|45| ). the center cavities between 
the two frequency-tunable cavities form a supercavity B27I1 . 
We note that our approach is also valid for the case A\ + A^. 




FIG. 10: (Color online) Circuit model of a frequency-tunable trans- 
mission line resonator integrated with a SQUID |31], 



IV. PHYSICAL REALIZATION OF 
FREQUENCY-TUNABLE SUPERCONDUCTING 
TRANSMISSION LINE RESONATORS 

In this section, we study several physical realizations of 
a frequency-tunable coupled-cavity array by using supercon- 
ducting transmission line resonators. In recent years, there 
have been several theoretical proposals and experiments on 
how to realize a frequency-tunable transmission line resonator 
(e.g., Refs. ll-Si). 

Typically, there are two physical mechanisms to tune the 
resonant frequency of a superconducting transmission line 
resonator. One method is to change the boundary condition of 
the electromagnetic wave in a transmission line. By changing 
the boundary condition, the effective wavelengths (also effec- 
tive frequencies) of the resonant modes are changed ll3ll - [33ll . 

Another method is to construct a transmission line res- 
onator by using a series of magnetic-flux-biased SQUIDs. 
Since the effective inductor of a magnetic-flux-biased SQUID 
can be tuned by changing the biased magnetic flux [34-36], 
the inductance per unit length of the SQUID array is control- 
lable. Therefore, the resonant frequencies of the modes in the 
SQUID array can be tuned by controlling the biased magnetic 
flux threading through the SQUIDs. 

Below, we present a brief review of these two methods al- 
ready used to obtain frequency-tunable transmission line res- 
onator. The original derivations of the two methods have been 
given in Refs. 113 ll 13511 . but for the sake of completeness of 
this paper, here we briefly review the main aspects of these. 



A. Ttining the frequency of a superconducting transmission 
line resonator: changing the boundary condition 

We briefly summarize the mechanism for frequency tun- 
ability of a superconducting transmission line resonator by 
controlling its boundary condition [31]. The lumped element 
circuit of a superconducting transmission line resonator with a 
symmetric SQUID, which is equivalent to a chain of identical 
LC circuits, is shown in Fig. [10] Here, <pj is the phase vari- 
able of the y'th node; C and L are respectively the capacitance 
and inductance of each LC circuit; <j> s \ and <p S 2 are the phase 



variables across the left and right Josephson junctions in 
SQUID, respectively; C s is the capacitance of one junctioi 
the SQUID. The SQUID is equivalent to a junction with 
effective Josephson energy E j(f) -2Ej cos (f/2). Here E 
the Josephson energy of one junction, / = 2n<$> x /<f>o, wh 
O x is the flux through the loop of the SQUID and Oo is 
magnetic flux quanta. <p s = (<p s i + <p S 2)/2 is the net ph 
across the SQUID. Note that here the self-inductance of 
superconducting loop is neglected. When ^ « 1, and 
charging energy and Josephson energy satisfy the condit 
E j(f) » C s [<&ol{2n)] 2 , then the SQUID can be approxima 
as a harmonic oscillator 13111 . For the phase variable of 
transmission line resonator, the wave equation reads ll3~lll 



ip (x, t) - v 2 (p" (x, t) = 0, 



(' 



where v = 1 / ^/CqLq. Here Co and Lq are, respectively, 
capacitance and inductance per unit length of the transn 
sion line resonator; and cp" refers to the second-order spa 
derivative. Note that here we have used the continuous v: 
able x instead of the discrete variable j. 

In terms of the relation between the electric current and 
phase I(x,t) — -^>o<p'(x,t)/(27iLo), the boundary condit 
for this system is IBlll 



/(0,f) = 0, 0Oo,O = «MO, 



( 



where /o is the length of the transmission line 1 
onator. The wave equation (l46b with the boundary cor 
tion ( |47| ) can be solved by assuming the solution <p (x, t) 
[A\ cos (kvt) + A2 sin (kvt)] cos (kx). The Euler-Lagrange 
equation for the phase variable <p s leads to the following dis- 
persion equation lf3lll 



2nY 
<6<J 



kl tan (JUo) = |-^| Ej(f) L cav - ^k 2 l 2 , 



2C S 
r 



(48) 



where L cav = LqIq and C cav = CqIq- The wave vectors k of 
the resonant modes in the transmission line resonator are the 
solutions of the dispersion equation ( l48l l. Since the effective 
Josephson energy E j (f) of the SQUID is tunable through the 
bias magnetic flux cD x , the wave vectors k can be tuned contin- 
uously by controlling <t> v . This approach l3lll to tune cavities 
could be used to tune the frequency of either one cavity or two 
cavities in our proposal. 



4 



(b) 



(a) 



^eff 



/JjfTL 

(C) 



FIG. 1 1 : (Color online) (a) Circuit model of a series array of 
SQUIDs. (b) Circuit model of a SQUID, which is equivalent to an 
effective tunable inductor in (c). 



where I c is the critical current of a single Josephson junc- 
tion, (pi (j = 1,2) are the phases across the two Joseph- 
son junctions. Introducing new phase variables <p — (<p\ + 
(pz)/2 and (pi - <pi = <p.x, then 7 = 7 C (<p x ) sin <p, where 
7 C (<Px) — 2I C cos(0 v /2). The phase can be expressed as 
<p — arcsin [7/7 c ((p x )]. An effective inductance L e ff can be de- 
fined [35] by the phase (p and current 7, 



^eff 



2tt7' 



which can be expressed as [35J 



7-eff (7, <p x ) = arcsin 

2tt7 \I c ((px) 



(50) 



(51) 



Obviously, the effective inductance of the symmetric SQUID 
can be controlled through two externally controllable param- 
eters: the biasing current 7 and the external biasing flux 
<t> v = Q>()4> x I(2ji). Thus the resonant frequency of the modes 
in a SQUID array becomes tunable because the inductance per 
unit length of the center conductor of the transmission line res- 
onator is controllable. This approach l35ll could also be used 
to tune the frequency of either one cavity or two cavities in 
our proposed system. This would allow the exploration of the 
effect predicted here. 



B. Tuning the frequency of a superconducting transmission 
line resonator: changing the effective inductance of a SQUID 

Following Ref. 135], Figure fTTT a) shows a device where 
the center of the resonator is composed of a series array of 
SQUIDs. A symmetric SQUID (in Fig.fTTrbT) is equivalent to 
an effective tunable inductance (shown in Fig. [TTT c)). From 
Kirchhoff's current law and the Josephson current-phase rela- 
tion Q, then 

7 = I c sin <pi + I c sin <p 2 , (49) 



C. Experimental implementation of our proposal 

Let us now provide some remarks on the experimental im- 
plementation of our proposal. In our model, the key elements 
are the frequency-tunable cavities, which have recently been 
realized experimentally. For example, in Ref. IB2I1 . the reso- 
nant frequency a> c of a transmission line resonator was tuned 
from 2jt x 4 GHz to 2k x 4.8 GHz (i.e., < \Aoj c \ <2nx 800 
MHz). If we choose a> c ~ 2n x 4 GHz, then < A < 0.2, if we 
choose a> c ~ 2n x 4.8 GHz, then -0.2 < A < 0. Similarly for 
Aq. In principle, the hopping coupling J between two nearest- 
neighbor transmission line resonators can be tuned 11611 . In 



10 



recent experiments (e.g., in Ref. l4lll ). the magnitude of the 
hopping interaction is / « 2n x 44 MHz w O.Olw,,. This hop- 
ping coupling 7 can be increased by using larger capacitors 
to connect two transmission line resonators. Therefore, this 
study seems to be within the reach of current (or near future) 
experiments. 

Compared to the method using a two-level atom as a con- 
troller, in Refs. Il4l 12711 . the present proposal avoids photon 
dissipation due to the spontaneous emission of the atom. 

It should be pointed out that we have neglected the change 
of the hopping coupling / between the frequency-tunable cav- 
ity and its nearest-neighbor cavities when the frequency of the 
FTC is tuned. In practice, this dependence exists. 

V. SUMMARY 

In conclusion, we have studied controllable single-photon 
transport and single-photon states in a one-dimensional 
coupled-cavity array (CCA) with one or two frequency- 
tunable cavities (FTCs). We found that, by adjusting the fre- 
quency of the frequency-tunable cavities, the coherent trans- 
port of a single photon in the CCA can be realized. We have 
also shown that there exist bound states in the CCA. 

For a CCA with one FTC, when the frequency of the FTC 
is larger than those of other cavities, there exists a bound state 
above the energy band of the CCA. When the frequency of 
the FTC is smaller than those of other cavities, there exists 



a bound state below the energy band. In these two cases, 
the bound states have even parity. Once the frequency of the 
FTC is given, the CCA has only one bound state. This re- 
sult is different from that of a CCA coupled with a two-level 
atom [14j,|27J], in which there exists two bound states at the 
same time, one above the energy band and the other below the 
energy band. 

For a CCA with two FTCs, in the two cases that the fre- 
quencies of the two FTCs are larger or lower than those of 
other cavities, there exist two bound states, one of odd parity 
and the other one of even parity. When the frequency detuning 
AqoJc of the two FTCs is very larger than the hopping coupling 
J between two nearest-neighbor cavities, there exist two kinds 
of resonant modes, one of odd parity, and another one of even 
parity. 
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